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The superfluid Reynolds number and the transition from potential flow to turbulence 

in superfluid ^He at mK temperatures 
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This comment is on Phys.Rev.Lett. 144, 155302 (2015) by M.T. Reeves, T.P. Billam, B.P. Ander¬ 
son, and A.S. Bradley “Identifying a superfiuid Reynolds number via dynamical similarity” where 
a new superfluid Reynolds number is introduced. This definition is shown to be useful in the data 
analysis of the finite lifetime of turbulence observed with an oscillating sphere in superfiuid helium 
at mK temperatures in a small velocity interval An = (v — Vc) just above the critical velocity Vc- 
The very rapid increase of the lifetime with increasing superfiuid Reynolds number is compared with 
the “supertransient” turbulence observed in classical pipe flow. 
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1. Introduction. In classical hydrodynamics the con¬ 
cept of the Reynolds number vD/v (where v is the flow 
velocity, D is the size of the object in the velocity field, 
and V is the kinematic viscosity of the liquid) is very 
useful because it reflects the similarity properties of the 
Navier-Stokes equation. For example, the drag coeffi¬ 
cient of a sphere moving through the liquid is a universal 
function of the Reynolds number only, irrespective of in¬ 
dividual values of v,D, and :/.[T] 

In a superfiuid the situation is more complex because 
there are three different velocities, namely the velocity of 
the superfiuid component Vs, the velocity of the normal 
component u„, and the motion of the quantized vortices 
with respect to Vg and Vn, the latter one was shown to 
be an intrinsic property being independent of the flow 
velocities, but only dependent on the interaction between 
both components. In superfiuid ^He those parameters 
lead to three different Reynolds numbers. [2] 

In superfiuid ^He below 0.5 K the normal component 
is almost absent, only a dilute gas of ballistically prop¬ 
agating phonons is left. Therefore, a hydrodynamic ve¬ 
locity Vn does no longer exist. In addition, it was shown 
in Ref. [2] that at these low temperatures the motion of 
the vortices in ^He is undamped and, hence, they are 
subject to turbulent motion. As a result, only one “su¬ 
perfiuid Reynolds number” Reg remains to be defined, 
namely by replacing the kinematic viscosity v by the cir¬ 
culation quantum of the vortices n = h/m (where h is 
Planck’s constant and m is the mass of a helium atom) 
because k has the same dimension m^/s as the viscos¬ 
ity: Rcg = vD/k.^ It describes the ratio of an over¬ 
all circulation to that of one vortex. In fact, in several 
theoretical investigations turbulence has been found for 
Reg » l.|l] 

In the new work by Reeves et al. [S] the two- 
dimensional Gross-Pitaevskii equation is investigated nu¬ 
merically in the vicinity of the critical velocity Vc for the 
onset of quantum turbulence. The central result of that 
work is the observation of dynamic similarity in the wake 
of a cylindrical object and its breakdown due to vortex 


shedding at a superfiuid Reynolds number defined as 


(The authors define here k = h/m; in the following we 
continue to use h/m.) Turbulence is observed in the 
wake at Rtg « 0.1, irrespective of cylinder size. 

2. The critical velocity Vc- In order to prove that 
(1) is a useful definition for the onset of turbulence in 
superfiuid helium, at first a determination of the critical 
velocity Vc is necessary. It is known that Vc depends on 
the experimental setup. 

Firstly, in experiments with steady flow the Feynman 
critical velocity Vc n/D is expected for vD/k ~ 1. 
Inserting this into (1) gives 



Only the zero point is shifted between both superfiuid 
Reynolds numbers. There is no difference when Reg ^1 
or u >> Vc- 

Secondly, in practice experiments with oscillatory flow 
are more easily performed, are more reliable, and in 
general yield sharp critical velocities. In the experiments 
with oscillating spheres [6], wires [7], and tuning forks 
[8] at various frequencies ui/2 'k ranging from ca. 100 Hz 
up to 30 kHz a critical velocity has been measured that 
scales as Vc ~ ^/kuj (neglecting at the moment numerical 
prefactors of order 1). Considering the oscillation ampli¬ 
tude v/lo as the relevant length scale D mvD/K^ 1 we 
obtain the experimental scaling. A rigorous derivation, 
however, is not available because there is still no theory 
of the transition to turbulence in oscillatory superflows. 

3. Transition to turbulence around an oscillat¬ 
ing sphere at Vc- Equation (1) is obviously relevant 
when the velocity amplitude is close to Vc, in particular 
when Av/vc 1. This is the case in our experiments 
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FIG. 1: (Color online) The normalized lifetimes r* = r/ro as 
a function of Av = v — Vc for the 119 Hz oscillator at 301 mK 
(blue squares) and the 160 Hz oscillator at 30 mK with 0.05% 
^He (black dots). Note the rapid increase of r* by 3 orders of 
magnitude over the small velocity interval of ca. 0.7 mm/s. 
The frequency, the temperature, and the ®He concentration 
have no effect on the data. The dashed curve is calculated 
from Eq.(3). 

with an oscillating sphere (radius R = 0.12 mm, frequen¬ 
cies 119 Hz and 160 Hz, Vc ~ 20 mm/s) at helium temper¬ 
atures below 0.5 K.[5] In a small interval of Av = {v — Vc) 
of few percent above Vc we observe that the flow pat¬ 
tern is unstable, switching intermittently between turbu¬ 
lent phases and potential flow. These patterns are easily 
identified because the drag force on the sphere is much 
larger in the turbulent regime than during potential flow. 
Recording time series at constant temperature and driv¬ 
ing force we analyze the distribution of the lifetimes of 
both phases. The turbulent phases have an exponential 
distribution exp(—t/r) and the mean lifetimes r increase 
very fast as we increase the velocity amplitude, namely 
as 

r = To exp (Au/ui)^, (3) 

where vi = 0.48 k/R, and tq = 0.5 s at 119 Hz and 0.25 
s at 160 Hz [9]. The velocity vi is of the same order as 
the Feynman critical velocity. The numerical prefactor of 
0.48 is due to various numerical factors which can safely 
be considered to be independent of the size of the sphere. 
Only the prefactor in Vc = 2.8 i/kw has not been proven 
to be independent because the size of the sphere was not 
varied. However, with oscillating wires nni the prefactor 
of Vc is 2.1, and this can be understood simply when 
comparing the maximum velocity at the equator of the 


sphere that is 1.5 times the fluid velocity far away with 
the maximum velocity at the circumference of a cylinder 
that is twice the velocity far away, therefore 2/1.5 = 1.33 
and 2.8/2.1 = 1.33. Hence, we expect the prefactor to 
have no significant dependence on the size of the sphere, 
if any. 

In Fig.l we plot the normalized lifetime r* = t/tq 
vs. Av. The salient feature is that r* is independent 
of the oscillation frequency, of the temperature, and 
is not affected by ^He impurities. The only frequency 
dependence is in tq. The exponential dependence of r on 
{Av/vi)‘^ may be described by Rice’s formula [TT1|T2] for 
Gaussian fluctuations crossing a given level. It is applied 
here to fluctuations of v crossing Vc where turbulence 
breaks down. The average number per unit time of 
crossings of a given level C below a mean at zero with 
negative slope is proportional to exp(—C'^/2cr^) with 
being the variance of the fluctuations. Postulating 
that the lifetime of a turbulent phase is ending at the 
first crossing of v of level Vc , we obtain (3) [T3]. Tq is 
interpreted to be twice the average time between two 
crossings of the level Vc [H], i.e., for tq = 0.5 s or 0.25 
s there are 4 or 8 crossings per second which means 
that the velocity of the superfluid around the oscillating 
sphere fluctuates at Vc with an average frequency of 2 Hz 
or 4 Hz, respectively. This could not be detected with our 
electronics. So far, we have no explanation for this result. 

4. Supertransient turbulence in '^He. Turbulence 
having a finite lifetime before breakdown to laminar flow 
is called “transient” turbulence, and in case of lifetimes 
growing faster than exponentially with Reynolds num¬ 
ber the term “supertransient” is used. m An example 
of supertransients was observed in classical pipe flow 
where the lifetimes of turbulent “puffs” have a double¬ 
exponential dependence on the Reynolds number. [T5] 

In our case we can replace Av by Rcg in (3) and set D 
— 2R to obtain 

r* = exp(ci?es)^, (4) 

with c = 1.04. In Fig.l the dependence of r* on 
Rcg is shown at the top axis (note that Rcg = 0 
corresponds to Av = 0 and r* = 1). It is obvious 
from Fig.l that our values of Rcs are all below 3, 
corresponding to maximum mean lifetimes r of about 
15 min and to Av/vc < 0.04. An extension towards 
larger velocities would require measuring times lasting 
much longer than our longest time series of 36 hours. 
Therefore, the validity of (4) cannot safely be extended 
to larger Rcg. Nevertheless, we can state that the fast 
increase of r* with Rcg is another example of supertran¬ 
sient behavior, here for the first time in a pure superfluid. 

5. Conclusion. The new superfluid Reynolds number 
is indeed a useful quantity for the investigation of the 
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transition to turbulence in oscillatory flow of superfluid 
helium just above the critical velocity Vc- The turbulent 
lifetimes grow very fast but without diverging at some 
particular superfluid Reynolds number. When the 
lifetimes begin to exceed the time of experimental 
observation, turbulence only appears to be stable, but 
actually it is supertransient chaos. 
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